Equations of Stellar Structure 


Stellar structure and evolution can be calculated via a series of 
differential equations involving mass, pressure, temperature, and 
density. For simplicity, we will assume spherical symmetry. 
First, let’s define the global variables 





r = the radial coordinate of the star 
M(r) 
L(r) 
P(r) 
p(r) 
X(r) = the mass fraction of the star that is hydrogen 
Y(r) 
Z(r) 
x;(r) 
u(r) 


7 = Total Mass of Star 

7 = Total Luminosity of Star 

R = the Radius of the star 

t = the independent time variable 





Obviously 
X+Y+2Z=1 


In addition, we will also define some local variables 


€ = the luminosity produced per unit mass of material 
Ky = the opacity (photon cross section) per unit mass of material. 


Conservation of Mass 


For a spherically symmetric star, the mass interior to some radius 
Ris 


R 
M(r) = f Arr? p(r) dr (2.1.1) 
0 
Written in terms of the differential, this is 
d 
M = 4rr° p(r) (2:12) 


However, over its lifetime, a star’s radius will change by many or- 
ders of magnitude, while its mass will remain relatively constant. 
Moreover, the amount of nuclear reactions occuring inside a star 
depends on the density and temperature, not where it is in the 
star. A better and more natural way to treat stellar structure 
is therefore to use mass as the independent parameter, rather 
than radius. Thus 


ce eee ee 
dM  4rr?p 


This is the Lagrangian form of the equation (rather than the 
Eulerian form). All the equations of stellar structure will be 
expressed in the Lagrangian form, and most of the parameters 
will be expressed in per unit mass, rather than per unit size or 
volume. 


(2.1.3) 


Conservation of Momentum 


If a star is in Hydrostatic Equilibrium, then the force of gravity 
at any point must be balanced by the pressure. Consider a shell 
at radius r: the mass per unit area in the shell is pdr, and the 
weight (per unit area) is —godr. This weight must cancel the 
pressure difference felt from one side of the shell to the other, 
Pinner ~ Pouter = IÈ dr. Thus 


dP 


can gree 2.2.1 
a gp (2.2.1) 


Since the force of gravity, g, is just GM /r?, then 





dP GM(r) 
= =P (2.2.2) 


So if we multiply by the conservation of mass equation, we get 
the Lagrangian form of hydrostatic equilibrium 











dP dr GM 1 
dp dM ( r? p) ae) ae 
or 
dP GM 
M Aas 28) 


Note: not all stars are in hydrostatic equilibrium. Specifically, 
some stars pulsate. For these stars, the pressure may not equal 
the force of gravity, and the shell will undergo an acceleration. 


Conservation of Energy 


Consider the net energy per second passing outward through 
a shell at radius r. If no energy is created in the shell, then 
the amount of energy in equals the amount of energy out, and 
d£L/dr = 0. However, if additional energy is released or absorbed 
within the shell, then d£/dr will be non-zero. Let’s define € as 
the energy released per second by a unit mass of matter. Then 





d 
= = 4rr°pe (2.3.1) 
or, in the Lagrangian form 
dL 


Note that e has three components. 

1) en, the total energy created by nuclear reactions 

2) €,, the energy directed into neutrinos, and 

3) €g, the energy produced or lost by gravitational expansion 
or contraction. 


Thus, 


dL 

dM 
In general, the contribution from nuclear reactions will always 
be positive, while the energy in neutrinos will always be lost 
from the system. (Supernovae are the exception to this rule.) 
The gravitational energy can be either positive or negative, and 
simply reflects the heat gained or lost in the shell by PdV work. 
For instance, our Sun will someday expand to ~ 100 times its 
size: this will require energy, and so each shell will absorb some 


En — €v + €g (2.3.3) 








energy and its eg will be negative. Similarly, a gravitational 


collapse will produce energy. Thus, 


dq 
Eg = ~ dt 


which, from thermodynamics (1.18), is 


dq dl’ odP 


oa Pt p dt 


(2.3.4) 


This can be written in a more useful form by recalling that 


Po 


Va~ 
: T pcp 





Thus 





oppi aE 
pg a ae P ae 


(1.19) 


(2.3.5) 


Thermal Structure 


The thermal structure of the star will depend on its method for 
transporting energy from its interior to its surface. For conve- 
nience, we will express the temperature stratification of the star 
in terms of the temperature gradient 


dln T 
= 241 
y (T) ( ) 


This temperature structure is then simply 





dT dP dT dP T 
na ae Ce) 


If the star is in approximate hydrostatic equilibrium, this simpli- 


fies to 
dT GMpT 


dr r2P 


or, in Lagrangian terms 


y (2.4.3) 





dT GMT 


dM‘ Anr4P ae) 


